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Mutual imitation games among artificial birds are studied. By 
employing a variety of mappings and game rules, the evolution 
to the edge between chaos and windows is universally confirmed. 
Some other general features are observed, including punctuated 
equilibria, and successive alternations of dominant species with 
temporal complexity. Diversity of species aided by the symboliza- 
tion of artificial birds' song are also shown. §1 Introduction The 

evolution of species in mutually interacting system has been a problem of 
interests in many branches of sciences. Is there trend to increase the com- 
plexity there? Is there any characterization for an evolutionarily favorable 
state? Determination of such a state is generally a complicated question, 
involving a certain balance among several factors, such as costs to adopt a 
strategy and behavior and gains by adopting it. The recent "dogma" for the 
characterization may be "the edge of chaos' scenario. This concept have been 
emphasized in many contexts such as cellular automata, boolean networks 
or CML 1^, 1^, 1^. To the authors' knowledge, however, there was no clear 
example presenting the evidence for evolution to the edge of chaos in the 
exact sense of dynamical systems theory. In the recent communication, we 
have proposed a minimal model for the evolution to the edge of chaos, based 
on the dynamical systems theory [^. The model has been motivated by the 
observed complexity of bird songs: It is known that a bird with a complex 
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song ( with many repertoires made from combinations of simple phrases) is 
stronger in defending its territory, as is observed by Krebs with the help of 
loud-speaker experiments [^l- So far the reason why a complex song is 
evolved is unknown, although there is a hypothesis that a complex song may 
give an impression of crowded population, thus is efficient for the defense0. 
This hypothesis, however, has no experimental evidences. On the other hand, 
there are some reports of observations that birds try to imitate each other's 
song at the defense of territory 0. Combining the above two observations, 
the following hypothesis has been put forward for the explanation of such 
complexity of a bird song: Birds join a battle for defense of territory through 
mimicking each other's song. In other words, if a bird A can imitate bird 
B's song better, then A can intrude B's territory and can share feeds. A 
is advantageous in the battle for the survival. It may be expected that a 
complex song is difficult to be imitated. Hence we may expect the evolution 
to a complex song. In spite of the lack of a direct experimental support for 
this scenario, the concept of mutually imitating game itself deserves much 
attention, as a novel problem of game theory and as a general framework 
for the evolution to a complex behavior. As an abstract game, we adopt an 
artificial bird (or player) whose song is generated by a nonlinear map. Birds 
play imitation games with each other. We assume that the winner of the 
game is fitter in the survival, which is taken into account by the population 
dynamics as a growth rate or by a simple replacement of losers by winners. 
Mutational changes of parameters for the song dynamics are introduced in 
the course of replication (See for details the modeling in section 2). We have 
to emphasize two advantages in our modeling here. 

• Since our artificial birds generate their songs by nonlinear mapping, 
the songs have complexity with real numbers. In most game theory, 
strategies are chosen from discrete sets, while the strategies (songs) 
here are chosen from continuously many sets. 

• Since our parameters have a clear meaning in the nonlinear mapping, 
the meaning of the edge of chaos is quite clear in the exact sense of 
dynamical systems theory. 

Since all players are labeled by continuous parameters, no clear speciation is 
given in advance. Judging from the result of simulations, however, it seems 
that we are entitled to divide players into several species. With high mutation 
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rates, the coexistence of several species has been observed. With the use of 
much smaller mutation rates, we have explicitly shown the evolution to the 
edge of chaos, in the exact sense of nonlinear dynamics. In addition, the 
evolution occurs towards the borderline between a window and chaos, not 
to the first onset of chaos. This observation leads to the hypothesis of the 
evolution to the border between observable and invisible chaos. Some other 
features emerge through the evolution, such as punctuated equilibria and 
abrupt, complex alternations of dominant species. In the previous report 
we have adopted 

• the logistic map as a song generator 

• euclidean criterion for imitation 

• "infinite dimensional" lattice (all-to-all game) 

• (see the meanings of these terms in section 2). 

Here we report studies on several variants of the model, to test the univer- 
sality of the results obtained previously. For example, we put birds on a 
2D square lattice, allowing only for fights among the players on the nearest 
neighbor sites. In most of cases, the results in the previous paper are re- 
produced, supporting the universality of the concept of the "edge of chaos" , 
and the edge of windows. Furthermore, we find several exotic phenomena 
such as the diversity (or coexistence) of species assisted by discrete symbol- 
ization of birds' songs. This paper is organized as follows. In the next section, 
we present details of our modeling. Several factors defining our models are 
described respectively. In section 3 we survey the results of the model with 
"Euclidean criterion" , parts of which were also reported previously . These 
results provide basic notions and framework for later sections. Characteristic 
features in the imitation games are shown there. We will give some intuitive 
reasons why the system evolves to the edge between chaos and a window. 
Moreover, it will be shown that our "edge of chaos" scenario is stable against 
choices of the topology of players. In section 4, results with the choice of 
symbolized criterion, i.e., "LR criterion" are reported. Such symbolization 
of dynamics is shown to lead to the increase of fluctuation, resulting the 
diversity of species. In section 5, we revisit the problem of the priority of 
the edge between chaos and a window, by changing the birds' song generator 
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(to the tent or the circle map). In particular, we will confirm the priority 
by simulating a system with the coexistence of two distinct types of song 
mappings with and without windows (i.e., the logistic and tent maps). In 
section 6, we will discuss complex dynamics of the imitation games by allow- 
ing for different dynamics for the "song generation" and imitation processes. 
Summary and discussions are given in section 7. §2 modeling Our model 

consists of the three stages, the song generation process, 2-persons' imitation 
game and reproduction. We will discuss these three procedures separately. 
1: Song generation (Song dynamics): As a "song", we use a time series 
generated by a nonlinear map. In the previous report [^, the logistic map: 



is adopted for this purpose, since it is one of the most investigated maps as 
a standard model in nonlinear physics. A bird, say, the i-th bird, possesses 
its own parameter a{i). We regard time series {xo{i),xi{i), ■ ■ ■) generated by 
the map as the i-th bird's song. For most examples of simulations here, the 
logistic map is again adopted, while the other maps for a bird's song generator 
are also adopted to test the universality of our results. Here simulations with 
the use of the tent map or the circle map will be given in §5. 
2: 2 person imitation game: Three factors are involved in defining the 
2-persons' game; topology for players, imitation , and game processes (or 
criterion for winners). Here we give detailed explanation of each process 
separately. 2. a) topology for players In the previous report 0, players are 
assumed to fight against all other players. In the sense of statistical physics, 
players live on an infinite dimensional lattice. Any spatial information is 
not included in this modeling. Here we also investigate cases where play- 
ers are put on a two-dimensional square lattice with the periodic boundary 
condition. For most simulations here, we adopt the lattice of 30 x 30. 2.b) 
imitation process: choice of an initial value for an imitating song: Each bird 
player i has to choose an initial condition, so that the time series by its own 
dynamics can imitate the other player's song better. Here we use the fol- 
lowing imitation process for simplicity. For given transient time steps Tj^j, 
a player (mimicker) modifies its dynamics with a feedback from the other 
player (singer): 




(1) 




(2) 
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By this dynamics, the player 1 adjusts its value x„(l) by referring to the other 
player's value (x„(2)). Here e is a coupling parameter for imitation process. 
After repeating this imitation process for Tj^j steps, the player 1 uses its own 
dynamics x„+i(l) = /i(x„(l)). In other words, the above process is used as a 
choice of the initial condition for the imitation of the other player's dynamics. 
The coupling parameter e also varies by players. However, the distribution 
of this parameter seems "irrelevant" , as far as we judge from our numerical 
results. Thus we will skip the discussion on this parameter hereafter. 
2.c): Game: After the player 1 completes the above imitation process, 2 
players 1 and 2 are decoupled and generate songs by their own dynamics. 
Then we calculate a quantity D(l,2) measuring the distance between the 
imitating time series {x„(l)}, and the singing one {a:„(2)}. By changing 
the role of two players, -D(2, 1) is also measured. If D(l,2) is smaller than 
-D(2, 1) the player 1 imitates better the other's song and wins this 2-persons' 
game ( and vice versa). Definition of the distance D{i,j) gives a criterion 
for the winner and may be crucial in the game. Here we adopt either of 
the following two criterions. 2.C.1) Euclidean distance: As the measure, we 
choose D(l,2) = I]m=i l^^mll) — a;m(2)P over certain time steps T. This 
choice is just a normal Euchdean distance. ( Of course, the distance should 
be cordial for the circle map, since x is on a circle.) 2.C.2) LR-code: 

It is often useful to symbolize the time series to some discrete sets rather 
than to use the real number x. The time series of the logistic map, for 
example, can be traced by 2 symbols, represented by R (x > 0) and L 
(x < 0) ||10|- Using such symbolization, the distance between two trajectories 



is measured by the number of unmatched symbols over some time steps. In 
the above symbolization for the logistic map, the distance is thus defined by 

T 

Dil,2)= Y.ISixmil),Xm{2)) (3) 

m=l 

with the notation IS{a, b) = 1 if ab < and otherwise. 3 Reproduction 

Fitness of the players are given by the results of games. Offsprings are 
reproduced according to the results of fights. We adopt either of the following 
two for this process : 3. a) Replacement: This might be the simplest way; the 
parameters of losers are replaced by those of winners at every game. We 
adopt this replacement rule for players on a 2D lattice. 3.b) Score: 

After each game, the winner gets a point W, while a loser gets L {W > L) 
(Both get {W + L)/2 in the case of draw.) After iterating a large number 
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of games, the population distribution is updated in proportional to the score 
of players, with the further restriction that the total population be constant. 
At both of these reproducing stages, we include mutational errors to the 
parameters; The parameters a and e are changed to a + S, e + S', where 
variables S and S' are random numbers chosen from a suitable distribution ( 
we use a homogeneous random distribution over [— /x, n] or Lorcnzian distri- 
bution {P{d) = + (5//^)^)}). The latter choice is often useful, since 
the former choice inhibits a large jump of parameters and often the parame- 
ter values are trapped at intermediate values, while the latter can provide a 
wider possibility of "species" . 

These three processes define our model. Because of the numbers of vari- 
ations, we will sometimes specify a model by sets of rules such as [lattice, 
logistic, LR, Score, 5e — 4], which means that players are on a 2D lattice, 
and generate songs by the logistic map, that the distance between two songs 
is measured by LR symbolic sequence, and that the population is updated 
in proportional to player's score with the mutation rate /x = 5 x 10^'^. 

Throughout this paper time step is measured in the unit of the total 
number of players N. In the lattice case one step is defined as the time 
necessary for one game for each player, i.e., totally N games for N players. 
In the all-to-all game, it is defined as the time required that each player plays 
a game for all othe players ( i.e., totally A^(A^ — l)/2 games for all players). 
§3 Basic Results: Simulations with Euclidean Criterion Let us start 

our presentation by giving numerical results with the Euclidean criterion 
for distance. First, we present a result of the game with high mutation 
parameters. The model is defined by [lattice, logistic, euclid, replacement, 
?]. In figure[l], we present the distribution of the logistic parameter a after 
10 time steps. It clearly shows the coexistence of two species; one is around 
a — 0.75 and the other is around a — 1.75. The former corresponds to the 
bifurcation point (period 1 period 2) while the latter to the edge between 
the period 3 window and chaos. The snapshot configuration is depicted in 
Fig [2], where dark dots mean players with a ~ 0.75 while bright dots give 
those with a ~ 1.75. This example already shows the typical character of 
advantageous states for some bifurcation points ( in other words, "edge of 
something"). 

To see the fitness by players we switch to use a model with all-to-all games 
[infinite, logistic, euclid, score, ?], presented in the previous communication. 
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since it deserves a prototype of our models. The average score of players is 
plotted as a function of the bifurcation parameter a by taking a high mutation 
rate, to allow for the existence of players with a wide range of parameters 
a[Fig 3]. We can see many peaks, which correspond to the bifurcation point 
from period-2 to period-4, and from period-4 to period-8, and the edge of 
period-3 window, period-5 window, period-4 window, and so forth. This 
score landscape is rugged as is often discussed in spin-glass type models in 
biology P, H . Here this landscape is not implemented as a model itself (like 
energy function in spinglass type models P]), but is emergent through the 
evolution. Indeed this landscape depends on the population distribution at 
the moment. To see the advantage of "the edge of chaos" in the score, we 
have plotted the score of birds as a function of the Lyapunov exponent A [Fig 
4]. As is shown there, the score has a clear peak around A = 0. Indeed this 
peak is supplied at the edge of period-3 window and chaos (a ~ 1-75), and 
the edge of period-4 window and chaos (a ~ 1.94). 

A high mutation rate, however, leads to a large fluctuation, which makes 
it difficult to identify each "species" precisely. Let us therefore concentrate 
on low mutation rate cases for the rest part of the present paper. In Fig. [5], 
temporal evolution of the average of the parameter a over all players is plot- 
ted. Plateaus are observed successively, providing an explicit example for 
punctuated equilibrium At the temporal domain with a plateau, the 



deviation from the average value is typically very small. Each plateau corre- 
sponds to the bifurcation points and the edges of chaos mentioned at the score 
landscape. Finally, the birds reach the parameter for the edge of period-4 
window (a ~ 1-94), and stay there [|. 

Thus we have seen the priority of "the edge of chaos" , and the evolution 
to it. To be more precise, we have to note that our edge state lies between a 
window and chaos. At windows, the logistic map can provide chaotic tran- 
sients before the dynamics settles down to a stable cycle. Thus the existence 
of transient chaos should be useful to imitate a dynamics of different nature. 
A window at a higher nonlinearity regime includes a variety of unstable cy- 



cles, as coded by Sharkovskii's ordering |]rUl. Therefore it can provide a larger 



variety of dynamics, as transients. This might be the reason why the edge 



^ The averaged score plotted as a function of the Lyapunov exponent A has a broad 
peak at A = 0, which extends to the region A < rather broadly. The score has a sharp 
drop at the side of A > 0, on the other hand. With the increase of T, the peak slightly 
shifts to A > while drops at A < get sharper. 
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of windows is strong in our imitation game. The above speculation suggests 
the importance of transient chaos, besides the edge of it, for the adaptation 
to a wide range of external dynamics. Let us incorporate again the lattice 
structure, by taking the model [lattice,logistic, euclid, score, ?]. Temporal 
evolution of the average logistic parameter value a is plotted in Fig[6]. Again 
we see the evolution to the edge of chaos (windows) after punctuated equi- 
libria. Thus the previous model with all-to-all game is worth of a good mean 
field theory to our problem. The lattice structure enhances the fluctuation as 
is seen by comparing the flatness of plateaus in Figs. (5) and (6) or noticing 
the lack of some plateaus here. To get some spatial information, we measure 
the distribution of the maximum connected cluster sizes of a species, since ad- 
missible parameter regions are separated enough to be regarded as distincted 
species. We have plotted the temporal changes of the maximal cluster sizes 
around the time steps for abrupt changes of the mean logistic parameter [Fig 
7] . Successive transitions of dominant species are clearly shown. An interest- 
ing question here is if the lattice structure helps the coexistence of species or 
not. In [Fig 7], both species with a ~ 1.75 and a ~ 1.86 increase their popu- 
lations together, which seems to suggest a kind of symbiosis. This example, 
however, is not so decisive to claim the spatial differentiation supporting the 
coexistence of species. Further studies are necessary. The examples in the 
present section show that the "edge of chaos" concept is valid irrespective of 
the topology of players. We stress again the significant role of the windows, 
which will be further addressed in section 5. §4 Symbolization induced 

complexity and diversity; LR criterion: Next we study the case with 
the use of the distance in terms of symbolic codes. The motivation for this 
choice is the examination of the effect of symbolization in communication 
codes: Does the symbolization lead to increase or decrease of diversity and 
complexity? The temporary evolution of the average logistic parameter is 
given in Fig[8]. In this case, advantageous players have such parameters that 
the symbolic sequence changes its pattern there, e.g., a = 1,1.3, It is 
not difficult to see that the dominant species lie again at window values in 
long intervals. Clearly, this example shows the stabihty of "edge between 
chaos and window" scenario against the choice of criterions in the 2-persons' 
game. Judging from the amplitudes in Fig [7] and [8], we conclude that the 
fluctuation here is further enhanced from the Euclidean criterion case. In the 
lattice version, the average values of a go up and down temporally with some 
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transient "disordered" intervals. No "final" parameter ( such as a ^ 1.94 
which is robust and is reached from any distributions in the Euclid criterion) 
exists here. The parameters switch (almost) forever. After staying at the 
edge between some window and chaos, the parameter switches to the edge of 
another bifurcation point or to an edge. In the example of Fig[8], successive 
transitions between a k, 1.75 (period 3) and a ~ 1.94 (period 4) are observed. 
Such complex alternations between "ordered states" are typically seen in the 
chaotic itinerancy [0, although the present switching process may be of 
stochastic nature (not of deterministic). Moreover, coexistence of different 
species (parameters) in the LR case is strongly enhanced compared with the 
Euclidean case with the same mutation rate. See Fig. [9] for snapshots of the 
distribution of players with respect to distinct sets of the logistic parameter 
a. The diversity is enhanced in the symbolic criterion. The increase of diver- 
sity is also seen in simulations with the all-to-all game. Roughly speaking, 
players with a ^ 1.75 (period-3 window) are rather robust and have occupied 
a large ratio of population. ( As the fluctuation is increased, a window with 
a larger interval (such as the period-3 one) is more robust, which is also true 
of a model with a high mutation rate with Euclidean criterion |^). Still, 
we have seen coexistence of other parameters such as a ~ 1.94 besides the 
above group. When the average is plotted with time, it fluctuates around 
a ~ 1.75 ~ 1.78, while the variance of the fluctuation remains large even 
in a very small mutation rate regime ( say 10~^). Summing up the present 
section, we have found that the symbolization enhances complexity of dy- 
namics and diversity of species. We have found coexistence of two groups, 
and successive switches of dominant species for ever. What is the origin of 
this diversity and complexity? One possible guess is as follows. By the binary 
symbolization, players with a certain range of the logistic parameter come 
to possess an identical temporal pattern in singing, whereas fine structures 
within the parameter interval are distinguished in the Euclidean (or ana- 
logue) measure. Thus a "strong" pattern of time series which was observed 
only in some very narrow window regime may be stabilized and exist in a 
wider parameter regime. The species, for example, at a ~ 1.63, can have a 
larger chance for the survival. By the above mechanism of symbolization, 
effective strength of players is further averaged out. Difference between two 
songs in the game is much weaker than the Euclidean one 0. Relationships 
between two players with different edges are more subtle, and no robust pa- 
rameter exists. The switching dynamics and the coexistence may reflect this 
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averaging effect by the symbolization. If tliis scenario of diversity is uni- 
versal, it suggests tlie importance of symbolization or discretization for tlie 
diversity. Why is the human language so diverse? Is the diversity caused by 
our symbolization ability from analogue vocal signals? §5 Universality of 

the priority of the edge between window and chaos Here wc examine 
the universality of the scenario in §3, by choosing a few different mappings 
for the song dynamics. In this section, we will simulate using the tent map 
and the circle map as the song generator. The results support our hypothesis 
on the priority of the windows, discussed in §3. 

Tent Map We choose the following tent map parameterized by one parameter, 

1 1 

Xn+l^a{--\--Xn\), (4) 

as the song dynamics. Thus each bird possesses two parameters a(i), e(i) as 
in the case of the game with the logistic map. In the tent map, there is no 
window structure. With the increase of a, chaos appears at a = 1, which is 
the only edge of chaos. For a > 1, no bifurcation structure exists. Numerical 
results of the evolution of the imitation game show that the average value 
of the parameter a evolves to 1, the onset of chaos for the tent map [Fig 
10, 11], irrespective of criterions. Since a = 1 is the only edge, this result is 
expected. For the logistic map, we have given a plausible explanation to the 
priority of windows in section 3. Is this edge of window stronger than the 
above onset of chaos in the tent map? Since no window exists in the tent 
map, it is interesting to study a system where birds with the logistic map 
and with the tent map can coexist and compete. Here we slightly modify the 
logistic map to 

^n+l = 2aXn{l — Xn) (5) 

SO that the variables (parameter a) in both maps take values in a same range 
[0,1] (a[0,2]). We put players on a 2D square lattice, and adopt the symbolic 
criterion (L for x < .5 and R for x > .5). Initially, a bird's song takes 

^ Another way of weakening the difference (but keeping ananlogue criterion) is the use 
of modified definition for distance measure; _D(1, 2) = km(l) — a;m(2)|^ with y <2. 

We have studied the case y = 1, which also enhances the fluctuation and instability of 
dominance of a species, similarly with the LR criterion. Such modification may be rather 
superficial, however, as is immediately seen by considering the extreme limit y ^ 0. 
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either logistic or tent map randomly with equal weights. Starting from the 
population with small a values, the mean value of a for logistic birds evolves 
to a punctuated equilibrium value a ~ 1 after few steps. (See Fig. [12]). 
At this equilibrium state, the ratio of the birds with the logistic map to 
those with the tent maps is kept almost constant ~ 8 : 1. Then the logistic 
parameter shows the abrupt change to the value corresponding to the period- 
3 window (or the period-4 window, depending on system size). During the 
change, the birds with tent maps are completely terminated, which means 
that the onset of chaos in the tent map is defeated by the windows in the 
logistic map. This result may deserve a support to our hypothesis on the 
priority of windows' edge. Circle map We have also studied the case with 
the use of the circle map 

Xn+i = Xn + a sin(27rx„) + d modi (6) 



as the song generator [|T^. Now each bird possesses three parameters a(i), d{i), e{i). 
We determine the rule for 2 persons' game by measuring the cordial distance 
between two trajectories of maps as noted previously. Numerical results 
again show abrupt changes of parameters with successive punctuated equi- 
libria. Synchronized temporal motions between average values of a and d 
are clearly seen in Fig. [13]. The average of d seems to take an "equilib- 
rium" value either ~ or ~ 0.5, resulting a shift of the origin of variables 
{xn{i)}- On the other hand, average a lies again at window values (to be 
precise, those with (i ~ 0) a ~ —2, —1.5, 1.5 etc. Summarizing this section, 
the significance of the window structures in the imitation game, discussed 
in section 3, is further confirmed by changing the generator for birds' songs. 
§6 Dual dynamics Before closing this report, let us present an example 



stressing clearly an aspect of dynamical complexity in the imitation game. 
Let the logistic map be again the song generator. This time, however, we 
assign two parameter values for the logistic map to each player so that it can 
use a different value in imitating and in singing. This is not an irrelevant 
complication of the model. Rather, there is no a priori reason to believe that 

To see the role of d, we also have carried out some simulations fixing d = and 
with same choices for other parameters. In this case, the magnitude of average a seems 
monotonously increasing with punctuated equilibria at windows ^ —5, —7 ■ ■ ■. Thus d is 
not irrelevant, but plays a role in suppressing the magnitude of the temporal variation of 
a. 
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a bird should use an identical parameter for the two processes. For imitation 
games with one parameter, the restriction of using an identical parameter in 
singing and mimicking may be a cause to achieve an optimal value of either 
singing or mimicking, and lead to an edge. In the two parameter game it 
may be possible that each process evolves to its own optimal value. Then one 
might expect that the results of the two parameter games' would be cither of 
the following: I. the parameters in singing and imitating evolves to the same 
"edge of chaos" value; or II. the parameter in singing reaches a value for a 
fully chaotic state (~ 2.0) while the imitating value also evolves to a high 
nonlinear region. Surprisingly, the observed result supports neither of these. 
In the Euclid criterion, the singing parameter again stays around a ^ 1.94 
the edge of window of period 4, while the imitating parameter stays around 
some other windows' edge ( e.g., a ^ 1.75) or other bifurcation points ( e.g., 
a ~ 0.75, bifurcation point from period-1 to 2), depending on the muta- 
tion rate and the topology of the game. Results with symbolic criterion are 
much more complex. Let us take an example with [lattice, LR, 5e-5]. The 
singing parameter switches among a ~ 1.77, 1.94,2, in synchronization with 
the changes of imitating parameter a ~ 1.1, 1.3, 1.77(Fig.[14]). This switch 
continues for ever as far as we have checked. The above result may be inter- 
preted as follows: Optimal values for singing and imitating parameters are 
not unique, but there may be several local minima. Thus optimal values for 
the singing parameter depend strongly on the population distribution of the 
mimicking parameter at that time, and vice versa. A small disturbance is 
enough to change both the optimal values and the population distribution 
completely. Synchronized changes between the two parameters reflect such 
a subtle balance. As stressed in section 3, the fitness function is not given 
a priori but is emergent through the evolution. For one-parameter games, 
however, the landscape is rather stable once the system evolves to the "edge 
of chaos" . In the dual dynamics game, the landscape, besides its ruggedness, 
varies in a quite complex manner. Thus the model provides an illustra- 
tion of a novel problem: evolution of a system with a dynamically changing 
rugged landscape. Summing up the present section, the evolution to the edge 
between chaos and windows is still valid in dual dynamics but the popula- 
tion dynamics is more complicated with successive alternations of dominant 
species, and synchronized changes between the two parameters. This exam- 
ple illustrates that static (or "equilibrium" ) characterizations are not enough 
for the imitation games, but the understanding of dynamics is essential to 
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them. §7 Summary and Discussions In this report, we have examined 

the evolution of imitation games. It is shown that the evolution to the "edge 
of chaos" is universally observed irrespective of topology of players, criterions 
for winners in a fight, and the choice of the map as a song generator. Besides 
the edge of chaos, we have also given a support to the priority of the edge 
of windows. It is clearly demonstrated by the dominance of the window's 
edge of the logistic map over the onset of chaos in logistic and tent maps. 
Though birds with both maps can coexist with parameter values around the 
onset of chaos, the players with the tent map are completely terminated with 
the abrupt change of the logistic parameter to the pcriod-3 window's value. 
Thus we propose an additional scenario here; the evolution towards the 
edge between chaos and window. This edge corresponds to the border 
between observable and invisible chaos. At a window value, there is topolog- 
ical chaos in dynamics. Although the final attractor is a periodic cycle for 
almost all initial conditions with a probability measure 1, there are chaotic 
orbits from non- measurable initial conditions (on a Cantor set). In connec- 
tion with this topological chaos, there are chaotic transients before an orbit 
is attracted to the periodic cycle. The existence of topological chaos assures 
a variety of unstable periodic orbits. Thus a player with a window parame- 
ter can imitate a large variety of periodic orbits, showing its dominance over 
players with periodic dynamics without topological chaos (a < 1.4011 •• • in 
our logistic map). By the transient chaos, the dynamics has an ability to 
imitate chaotic time series roughly up to the length of transients, which di- 
verges at the edge. As a song generator, dynamics at a window parameter 
can provide a large variety of orbits as transients. Since the transient length 
diverges at the edge, a high variety of songs is maintained. Thus generated 
songs are not easily imitated by periodic or chaotic dynamics. Summing up, 
our window's edge scenario is based on the abihty of creation of complexity 
with topological and transient chaos. Transient chaos has a potentiality to 
adapt a wide range of external dynamics, while the orbit at the attractor 
is not complicated. The transient chaos may be important in a wide area 
of biological information processing, and our window's edge scenario may 
be applied universally in the evolution and adaptation. Besides the char- 
acterization of the evolutionary advantageous states, we have discussed the 
origin and maintenance of diversity of species, in a symbolized criterion and 
in dual dynamics for songs and mimicry. If there is a given fixed fitness 
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landscape, it would be reasonable to expect that the fittest species would 
dominate the world. Indeed, without any mutations or changes of environ- 
ment the diversity of species may not be sustained. There can be several 
possibilities in the origin of diversity. A simple answer to the origin is the 
spatial differentiation mechanism. Indeed our simulation shows coexistence 
of some species in a 2-dimensional topology. The diversity is increased with 
the use of spatial game. However, a more interesting discovery in the present 
paper is symbolization induced diversity. By adopting a symbolized cri- 
terion, the diversity is enhanced drastically. At some stage of information 
processing in the brain, symbolization is often adopted. Thus we may expect 
that the symbolization induced complexity can be one origin of diversity in 
signal and language. In the symbolized criterion, the population dynamics 
is also complicated. It shows successive punctuated equilibria for ever. The 
complexity is further enhanced with the choice of dual dynamics for songs 
and imitations. A variety of window's edge states appear successively, pro- 
viding temporal complexity. In a large system size, this complexity is easily 
expected to lead to spatial diversity, since the successive changes of domi- 
nant species cannot be synchronized over all lattice regions. The diversity 



induced by temporal chaos is recently discussed as homeochaos scenario ^ 
where the population dynamics of species shows weak and high-dimensional 
chaos. Our window's edge scenario shares the notion of "weak chaos" with 
the homeochaos scenario. In the homeochaos, we have also seen successive 



switches among ordered states, noted as chaotic itinerancy [Ig, |T6l . Detailed 
study of dynamical mechanisms of the successive switching in our case will 
be necessary in future. Our imitation game provides a universal route to the 
evolution to complexity. Evolution pressure of complexity to escape from 
imitation can be also conceptually used in many examples in the biological 
or social evolution, besides our original motivation to a bird song. Such ex- 
amples may include the evolution of a communication code ( "secret code"). 
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Did God tell us the symbolization to destroy the Tower of Babel? 
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Figure Caption Figure 1 Snapshot distribution of the logistic parameter of the 

players after 10 time steps. Mutation parameter /i is chosen to be 0.001, and 
lattice size is 30 x 30. Two peaks corresponding special points for the logistic 
map are observed (see text). Unless otherwise mentioned, we set Tj^j = 
T = 30 throughout the present paper. Figure 2 Snapshot pattern of the 

configuration of distribution of the logistic parameter, corresponding to Fig.l. 
Here a bright (dark) pixel corresponds to a player possessing a ~ 1.75(0.77), 
respectively. Figure 3 Emergent landscape: Average score for the players 

with parameters within [oj, Oj + A] is plotted for aj = — 1+ixA, with the 
bin size A = 0.001. We have adopted = 10, L = 1, Timi = 255, and T = 
32. Simulation is carried out with the rule [infinite, logistic, Euclid, score, 0.1]. 
starting from the initial parameter a = 0.6 and e = .1. Sampled for time steps 
from 1000 to 1500, over all players ( whose number is fixed at 200). (adapted 
from 0). Figure 4 Average score of the game vs. Lyapunov exponents. 

Simulation is carried out with the rule [infinite, logistic, Euclid, score, 0.02], 
Timi = 255 and T = 32. We have adopted W = 10, L = 1. starting from 
the initial parameter a = 0.9 and e = .1. Average scores are obtained from 
the histogram of Lyapunov exponents, for which we use a bin size of 0.01 
for — 1 < A < 1, while it is set at 0.1 for A < — 1 ( since the sample there 
is rather sparse). Sampled over time steps from 500 to 750 over all players 
( whose number is fixed at 200). Figure 5 Two examples of the temporal 

evolution the average parameter a: Simulation is carried out with the rule 
[infinite, logistic, Euclid, score] , and with mutation rates 5 x 10~^ and 8 x 10~^. 
Timi = 200 and T = 50. starting from the initial parameter a = 0.6 and 
e = .1. Average of the parameters a over all players are plotted with time. 
The total population is fixed at = 200. Figure 6 Temporal evolution of 

the average logistic parameter a. (as is written, 1 time step = 900 2 persons' 
game), [lattice, logistic, Euclid, replacement, = 8.7 x 10~^] Initial values 
for a are chosen randomly from [0.6, 1.1]. Punctuated equilibria at a ~ 0.77 
and 1.94 are clearly seen. The plateau at 1.94 continues after 3500 steps, 
as far as we have observed. Figure 7 Temporal evolution of 4 connected 

clusters in the vicinity of the change of the average logistic parameter in 
Fig.6. Around the time steps [2700, 2800], two clusters (~1.74 and ~ 1.86) 
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seem to grow simultaneously. Figure 8. Temporal evolution of the average 

logistic parameter a. [lattice, logistic, LR, replacement, mu = 1.25 x 10~^]. 
Alternations between the states at a ~ 1.77 and ~ 1.94 are seen. Figure 

9. Successive snapshot distribution of the logistic parameter at time steps 
1000,5000,7000. The parameters of the model arc taken identical to those 
in Fig. 8. Figure 10 Temporal evolution of the average tent parameter. The 
rule is given by [30 x 30 lattice, tent, EucHdean, replacement, 6.1 x 10~^]. 
Initial values for the tent parameter a are chosen randomly from [0.3,0.4]. 
The deviation from 1 is very small after establishing the equilibrium. Figure 

11 Temporal evolution of the average tent parameter with the LR criterion. 
The rule is given by [30x30 lattice, tent, LR, replacement, 8 x 10^^]. We have 
also simulated smaller fi cases, but the fluctuation remains still larger than 
that from the Euclidean criterion (given in the previous Figure). Figure 12 

Temporal evolution of the average value of logistic parameters a (thin line) 
and the population of adopting a logistic map song. Here the lattice size 
is 32 X 32 (totally 1024 players) and /j, = 7.5 x 10~^. Around the time 
step 4000, birds with logistic maps cover the total lattice and those with 
the tent maps are terminated. Figure 13 Temporal evolution of the average 

a (thin line) and ^(broken line) values. [30 x 30 lattice, circle, Euclidean, 

replacement, 5 x 10~^]. Figure 14 Temporal evolution of the averages of the 
logistic paremters for singing (broken) and imitating (thin) in dual dynamics, 
with LR criterion. The lattice size is chosen to be 30 x 30, and = 7.5x 10~^. 
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